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ABSTRACT 

Context. The physical conditions in massive dense cores (MDCs) forming massive stars and clusters, are not well 
constrained. Observations are lacking to confront the theories. Recently an extensive study has been started towards 
the most massive and youngest cores in the Cygnus-X molecular complex and as a first result showed exceptional 
fragmentation properties in a sample of five cores where individual massive protostars have been recognized. 
Aims. We study here the kinematic properties of dense gas surrounding massive protostars in these five cores in order 
to investigate whether turbulent support plays a major role in stabilizing the whole core against a rapid fragmentation 
into Jeans-mass objects. Alternatively, the observed kinematics could indicate a high level of dynamics suggesting that 
the cores are actually not in equilibrium and dynamical processes could be the main driver to build up the final stellar 
masses. 

Methods. We present IRAM 30m single-dish (H 13 CO + and HCO + ) data and IRAM Plateau de Bure Interferometer high 
angular-resolution observations of dense gas tracers (H 13 CO + and H 13 CN) to reveal the kinematics of molecular gas at 
scales from 0.03 to 0.1 pc. 

Results. Using radiative transfer modeling we show that H 13 CO + is depleted within the envelopes of massive protostars 
and traces the bulk of material surrounding the protostars rather than their inner envelopes. H 13 CN shows a better 
correspondence with the peak of the continuum emission, possibly due to abundance anomalies and specific chemistry 
in the close vicinity of massive protostars. Analyzing the line-widths we show that the observed line-dispersion of 
H 13 CO + at the scale of MDCs is smaller than expected from the quasi-static, turbulent-core model. At large-scales, 
global organized bulk motions are identified for 3 of the MDCs. At small-scales, several spectral components are 
identified in all MDCs showing filamentary structures and intrinsic velocity gradients towards the continuum peaks. 
The dynamics of these flows show diversity among the sample and we link this to the specific fragmentation properties of 
the MDCs. Altogether this points to different initial conditions towards CygX-N3 and -N63 as compared to CygX-N12, 
-N48 and -N53, which may represent different evolutionary stages. 

Conclusions. No clear evidence is found for a turbulence regulated, equilibrium scenario within the sample of MDCs. We 
propose a picture in which MDCs are not in equilibrium and their dynamics is governed by small-scale converging flows, 
which may initiate star- format ion via their shears. We suggest that dynamical processes are linked to the formation of 
proto-clusters and high-mass protostars. 
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1. Introduction 

Several models have been proposed to explain what pro- 
cesses dominate and form structures within giant molecu- 
lar complexes, thus determining the origin of star- forming 
cores and ultimately of stars. After early scenarios based 
on quasi- static evolution of molecular clouds, primary reg- 
ulated by the interplay between gravitation and magnetic 



mated by tne interplay between gravitation ana magnetic 
fields (e.g. |Mouschovias fc Spitzer| |l976l |Shu et al.|fl987| 



Mouschovias fc Ciolek|1999p , the importance of turbulence 
has been recently emphasized and taken into account (see 
Mac Low fc Klessen|2004l|McKee fc Ostriker|20"07|and refer- 



ences therein). For high-mass star formation, the consider- 
ation of the relative importance of these different supports 
against gravity has led to two main alternate views. In the 
first scenario molecular clouds are close to be in equilibrium 
and are evolving on long time- scales (of the order of a few 
dynamical crossing times, e. g. Larson||1981 McKee||1999|). 
Since turbulence decays within a dynamical time (e.g. Mac 



Low et al.| 1998} |Padoan fc Nordlund||i999| ), it needs to 
be replenished over such long time-scales (e.g. by stellar 
winds, supernoyae explosions, ou tflows or other feedbacks, 
|Norman fc Silk|[l980| |McKee|[l989[ ) in order to maintain a 
quasi-static state. As a result, turbulence can regulate the 
formation of structures within molecular clouds on large 
scales, while at small scales isotropic, supersonic turbulence 
could be the main support of cores a gainst collapse, thus 
complementing their thermal pressure ( |McKee fc Tan|2002| 
2003 ) . It naturally offers a description for the formation of 



high-mass stars as a scaled-up scenario of the low-mass, 
quasi-static star- formation process. Accretion rates of up 
to ~ 10 -3 M yr _1 are reached for these turbulent, massive 
cores, which is high enough to overcome radiation pressure 
and enables the protostar to collect a mass larger than 8- 
10 Mq. The turbulence regulated, quasi-static scenario has 
fundamental predictions: 1.) the existence of massive cores, 
which are not yet forming stars (an equivalent of low-mass 
pre-stellar cores); 2.) a high-level of isotropic, supersonic 
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turbulence on the size-scale of massive cores; 3.) the col- 
lapse of massive cores directly into a single (or close bi- 
nary) massive object. To be more precise, gravitationally 
unstable density fluctuations are predicted, but with a low 
number of gravitationally bound fragments within a single 



core. 



Alternatively, clouds and clumps may be unstable and 
fast evolving structures with high-level of dynamics (e.g 



Vazquez-Semadeni et al.|2002||Klessen fc Hennebelle|2QQ9[ ) . 
Briefly, large-scale turbulent flows create structures by 
shock-dissipation and, as shown by numerical simulations, 
supersonic turbule nce fragments the gas efficiently in very 



short time- scales (jPadoan et al. 2001 Vazquez-Semadeni 



et al.||2007 ). This picture of gravo -turbulent fragmentation 



naturally provides the seeds for star-formatio n by creating 
gravi t ationally bound density fluctuations ( |Klessen et al.| 
2005] |Dib et aL||2007| ). In hydrodynamic simulations this 
results in fragments with mass around the local Jeans mass 
(~1 M ), and thus, it is efficient to produce low mass- 
fragments, while there is a significant lack of massive ones 
flBonnell et aT1|2007| [Dib et al.||2010b[ ). Bonnell et al. pro- 
posed as a solution that seeds may continue to accrete ma- 
terial from regions which were originally not bound to the 
protostellar envelope. Thi s scenario is therefore introduced 
as competi tive accretion ( |Bonnell et al. 2001 Bonnell & 
Bate 2006). Massive seeds continue to grow in mass via two 
mechanism: benefitting from the channelling of material on 
the common cluster potential well or if the relative veloci- 
ties are too low, the accretion processes are dominated by 
the tidal field of the protostar itself. 



According to Hennebelle & Teyssier (2008) magnetic 
Gelds can lower the level of fragmentation resulting in only 
a few more massive objects, which continue to gain mass 
via channelling of material along the field lines. Note, how- 
ever, that strong perturbations (that must be seeded in 
the clumps/cores in order to allow them to fragment in 
the presence of magnetic fields) may also naturally be re- 
produced in the gravo-turbulent fragmentation of magnet- 



ically superc ritical and nearly critical cloud s as well (Dib 
et al.||2010a[). In the analytic description of Hennebelle h 



Chabrier| ( |2008p (magnetic field included) , the interplay be- 
tween gravothermal and gravoturbulent fragmentation can 
naturally reproduce the observed frequency of massive stars 
- depending on the Mach- number within molecular clouds. 
The fundamental predictions of the dynamical view are, 
that 1.) non-isotropic turbulence/velocity fields should be 
present with dynamical patterns; 2.) a high-level of frag- 
mentation is expected with the most massive stars forming 
in the center of clusters (since they are the most efficient 
in competing for mass). 

We intend to put observational constrains on the key pa- 
rameters of current star-formation theories by performing a 
systematic, high angular-res olution study of massiv e dense 
cores (MDCs) in Cygnus-X. Bontemps et al. ( 2010| ) (here- 
after Paper I) presented 1mm and 3mm continuum maps of 
6 IR-quiet MDCs, which were obtained by the Plateau de 
Bure Interferometer. Twenty three fragments were detected 
in a sample of six cores, which challenge the turbulence 
regulated, monolithic collapse view - according to which 
these cores should be much less fragmented. However, the 
most compact MDC, CygX-N63, stands out from the pic- 
ture since it contains only a single, compact object, and 
could therefore be a representative of a single massive pro- 
tostar with an envelope mass of ~60 M . In addition, a 
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Fig. 1. Spitzer /IRAC 8 fim map (Hora et al J2009| showing 



a part of Cygnus-X North. Black contours show integrated 
intensity map of the N2H + (J= l-0) line obtained with the 
FCRAO flSchneider et al.||2010| ). Red stars indicate the po- 
sition of the sample of MDCs. Green triangles ind icate the 
position of Hn regions ( Downes fc Rinehart||1966 ). 



total of 8 fragments are found to be precursors of OB stars 
with envelope masses ranging from 6-23 M . We found that 
the level of fragmentation within these MDCs is higher than 
expected in the quasi-static formation scenario, but lower 
than predicted by a pure gravo-turbulent scenario. Also, 
we compared the mass of the cores with the total mass 
contained in the fragments, which is 28, 44 and 100% for 
CygX-N12, CygX-N53 and CygX-N63, respectively. This 
efficiency of turning the cloud mass into protostars is sig- 
nificantly higher tha n what has been observed for clusters 
( |Lada fc Lada||2003|) and predicte d for low-mass stars and 
clusters QMatzner fe McKee||2000| ). 

The origin/physical reason of such diverse fragmen- 
tation properties is unknown, and can only be eluci- 
dated with kinematic studies using molecular lines. Here 
we present a study of high density molecular trac- 
ers using single-dish observations of HCO + (J=1-0) and 
H 13 CO + (J=1-0) line emission from a sample of MDCs to- 
gether width high angular-resolution, interferometric ob- 
servations of H 13 CO+(J=1-0) and H 13 CN (J=l-0) lines. 
Th ese molecules have a hi gh critical density (>10 4 cm -3 , 
e.g. |Bergin fc Ta falla 2007) and are known to be tracers of 
dense gas, therefore allowing to probe the dynamics of gas 
within MDCs. 



2. Source selection and observations 

2.1. I R- quiet MDCs in Cygnus-X and sample selection 

The source selection dates back to the 1.1mm MAMBO 



survey of Cygnus-X ( Motte et al.||2007 ), from which we se- 
lected the most massive IR-quiet MDCs. All cores are com- 
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pact with a similar size as nearby low-mass cores (~0.13 
pc), but they are 10 times denser (with an average density 
of 1.9 x 10 5 cm -3 ) and 20 times more massive (~ 60- 
200 M ). Therefore one would expect them to be actively 
forming stars, which is in fact confirmed by outflow emis- 
sion traced in SiO. However, the lack of strong mid-IR and 
free-free emission indicates that they must be in an early 
stage of their evolution. 

As shown in Figure [I] two of the selected cores (CygX- 
N48 and CygX-N53) are located in the DR21 filament 
which is the most massive and densest region of Cygnus-X, 
and is a well-known region of massive star formation (see 
e.g. Schneider et al. ( 2010[ ) and references therein). CygX- 
N48 corresponds to the submillimeter source DR21 (OH)-S, 
while CygX-N53 is situated close to W75-FIR3 flChandler 
et al.||1993[ ). The 3 other selected MDCs have been discov- 
ered by the MAMBO survey and are situated in more iso- 
lated, but still prominent molecular clumps bright in N 2 H + 
line emission, which equally traces cold, dense gas. CygX- 
N3 in the western part is located close to DR17 which is 
an Hn region most prob ably excited by two OB cluster s 
(clusters # 12 and 14 in |Le Duigou fc Knodlseder||2002 ) . 
These clusters are shaping the cloud and CygX-N3 corre- 
sponds to the tip of a pillar-like cloud. CygX-N12 is also 
located in a cometary shaped cloud, probably influenced by 
the same OB clusters. However the cloud seems to be less 



comp ressed from outside than CygX-N3 (Schneider et al. 



20Q6D - Finally, CygX-N63 is placed in the south of DR21 in 
the DR22-DR 23 filament. We adop t a distance of 1.7 kpc 



for Cygnus-X ( Schneider et al.|2 006) so that 1" corresponds 
to 1700 AU. 



2.2. Observations 

2.2.1. Observations with the IRAM 30m telescope 

OTF maps of spectral lines of CygX-N3, -N12 and -N63 
were observed on 5 June 2007 with the IRAIV0 30m tele- 
scope using the A100 receiver with the VESPA correlator 
for H 13 CO+ (J=l-0) at 89.188 GHz and HCO+ (J=l-0) 
at 86.75 GHz. The average system temperature was ~85 
K. Sources CygX-N48 and -N53 were observed between 2- 
5 June 2007 as part of a large-scale O TF mapping of the 
DR21 filament ( |Schneider et al. 2010), and have average 
system temperatures around 110 K. The data were cor- 
rected for the main beam efficiency of 0.78 and have a 
velocity resolution of 0.135 km s" 1 . The HPBW at this 
frequency is ~28". 



2.2.2. Observations with the Plateau de Bure Interferometer 

We used the IRAM Plateau de Bure Interferometer (here- 
after PdBI) to perform high angular-resolution observations 
at 1mm and 3mm. A more detailed description of the obser- 
vations can be found in Paper I| We obtained line emission 
measurements of H 13 CO + (J=l- 0), and H 13 CN (J=l- 0) 
at 86.75 and 86.34 GHz, respectively, with a velocity resolu- 
tion of 0.134 km s _1 . The observations were done in track- 
sharing mode with two targets per track. The D configura- 
tion tracks were obtained between June and October 2004 
with 5 antennas with baselines ranging from 24 m to 82 m. 



The C configuration tracks were obtained in November and 
December 2004 with 6 antennas in 6Cp configuration with 
baselines ranging from 48 m to 229 m. We mostly used as 
phase calibrator the bright nearby quasar 2013+370 and as 
flux calibrator the bright evolved star MWC349 which is 
located in Cygnus-X. 

2.3. Data reduction and cleaning 

We used the GILDAS softwar^E] for the data reduction 
and analysis of single-dish and interferometric data. The 
interferometric H 13 CO + dataset was reduced both with 
and without adding short-spacings information. The short- 
spacings information was obtained with the IRAM 30m 
telescope and implemented with the standard techniques 
of GILDAS. No zero-spacings information was available 
for the H 13 CN dataset. For the cleaning procedure we 
used a naturally weighted beam to search for clean com- 
ponents within a circle. The resulting parameters of the 
data-reduction are summarized in Table [I] All maps are 
corrected for primary beam attenuation. 

3. Results 

As a continuation of the interferometric mm-continuum 
studies of Paper I, we focus here on a study of the distri- 
bution and the kinematics of dense gas traced by HCO + , 
H 13 CO+ and H 13 CN in our sample of MDCs. We first 
present single-dish maps of HCO + and H 13 CO + , on which 
we base later our radiative transfer line modeling in Sect. 
4.1.2 Then, spectra from both single-dish and interfero- 
metric maps are shown, followed by line integrated maps of 
H 13 CO+ and H 13 CN. 



3.1. IRAM 30m H 13 CO+ and HCO+ observations 

Spectral line maps of the IRAM 30m observations of the 
five MDCs are shown in Figure [2] a). We obtained a map of 
line emission of ~1 arcmin 2 of the 3 isolated MDCs (CygX- 
N3, -N12, -N63). For Cyg X-N48 and -N53 we cut out maps 
from a larger mapping of Schneider et al. (2010) with the 
same dimensions. 



1 IRAM is supported by INSU/CNRS (France), MPG 
(Germany) and IGN (Spain). 



Both lines are well detected over the extent of MDCs 
at a size-scale of ~0.1 pc. The large-scale distribution of 
emission indicates a rich gas content in all cores. For the 
isolated cores (CygX-N3, -N12, -N63), the distribution of 
H 13 CO + emission is centrally peaked and confirms the core 
as a coherent entity. For the fields located in the DR21 fila- 
ment, the distribution of emission follows the main pattern 
of the filament and therefore shows a strong north-south 
elongation. 

The H 13 CO + line shape is a rather clear gaussian in 
CygX-N53 and -N63, while in the other fields it shows a 
nearly gaussian profile. Blending due to several components 
and broad line-wings are remarkably present in CygX- 
N3 and CygX-N12. Towards CygX-N48 a red-shoulder is 
present in the spectra, while the stronger component is cen- 
tered on the self-absorbed profiles of HCO + indicating two 
strong components in the spectra. 

The optically thick tracer, HCO + , shows blue-shifted 
asymmetry in four MDCs of our sample. In CygX-N3 this 
profile changes into a reversed inf all-profile, a red-shifted 



See http://www.iram.fr/IRAMFR/GILDAS/ 
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Fig. 2. a,) Spectral line maps of HCO + (black) and H 13 CO + (red) obtained with the IRAM 30m telescope and regridded to 
full-sampling. The velocity range is 10 to 20 km s _1 for CygX-N3, -N12, -15 to +5 km s -1 for CygX-N48 and CygX-N53, 
410 to km s _1 for CygX-N63. The intensity of the optically thin line (H 13 CO + ) is degraded by a factor of 2-4. Dashed 
squares indicate the size-scale of the right panel, b) Spectral line maps of H 13 CO + PdBI combined with IRAM 30m, 
where the spectra-cubes were regridded to 3" resolution. Velocity ranges are 10 to 20 km s _1 for field CygX-N3,- N12 
and -10 to km s _1 for CygX-N48, -N53 and -N63, respectively. The grey-scale background image is the 3mm continuum 
emission obtained si multaneou sly with the PdBI (from 3a to75 a). (The grey dashed contours show the same logarithmic 
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Table 1. Observation and data reduction parameters of the interferometric data. We give beam-sizes, position angle and 
noise levels of H 13 CO + for the combined dataset of PdBI+30m. For H 13 CN, parameters of the interferometric data are 
shown. 
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Fig. 2. - continued. 



peak appears over a part of the map. The dip is at the 
same velocity where the optically thin line peaks, which 
is commonly int erpreted as an exp anding envelope or an 
outflow feature (Myers et al. 1996). However, ~15" away 
from the center the profiles turn into a standard infall pro- 
file. Such blue-shifted asymmetric profiles are observed to- 



wards CygX-N48 and -N53. In all these cases the dip of 
self-absorption is located at the peak of the optically thin 
line indicating infall. Interestingly, CygX-N12 and espe- 
cially CygX-N63 show only a slight asymmetry in their 
spectra. Since the ratio of H 13 CO + and HCO + line inten- 
sities (7, 5.5, 5.5, 1.9, 3.7 for CygX-N3, -N12, -N48, -N53, - 
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Fig. 4. a) Integrated intensity maps of H 13 CO + (J=1-0) obtained with PdBI and combined with observations from the 
IRAM 30m telescope (the scaling goes from -2a to 35a for field N63 and N12, and from -2a to 60<j for all other fields). 
Contours show the 3mm continuum emission, where levels are the same as in Paper I Crosses indicate the location of 
1mm continuum sources reported in |Paper~l| Dashed circles indicate the primary beam of the PdBI. Grey ellipses in 
the left lower corner show the synthesized beam, b) Interferometric line integrated intensity maps of H 13 CO + (J=l-0). 
c) Integrated line intensity maps of H 13 CN (J=l-0) (over all hyperfine components). 
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Fig. 3. HFS line fitting of the H 13 CN spectra, which were 
averaged within 50% contours of the peak emission. Green 
line indicates the result of HFS line-fitting assuming LTE 
conditions. All hyperfine-components are detected towards 
all of the cores. 



N63, respectively) indicate that the emission of HCO + must 
be optically thick, we interpret this b lue-shifted asymme- 
try as infall signature. In Section [4. 1.2| we present radiative 
transfer modeling of these line profiles. 

3.2. IRAM PdBI H 13 CO+ and H 13 CN observations 

3.2.1. IRAM PdBI H 13 C0+ spectra maps 

Maps of H 13 CO + spectra obtained with the PdBI and pro- 
cessed with zero-spacings are presented in Figure [2 b). 
We regridded the spectra-cubes to a fully sampled syn- 
thetic HPBW of ~4". These maps of spectra taken at high 
angular-resolution resolve several distinct spectral compo- 
nents. For fields CygX-N3 and -N53 individual spectral 
components are found with a separation of 0.5-2 km s _1 
(4x-15x the spectral resolution). They appear clearly off 
the main position of the continuum peaks and can be sep- 
arated into two strong components each of them showing 
a shift in position and in velocity and tend to blend into 
a single line towards the continuum peaks. In fields CygX- 
N12 and -N48 the spectra are more complex, especially in 
CygX-N12 the spectra are flattened, which makes the sep- 
aration of the components difficult. Similarly in CygX-N48 
the two strong components seem to be blended with other 



components and in the central region broad lines are seen. 
(Note that for CygX- N48 a large- scale sub-f ilament at 
km s _1 is reported by |Schneider et aL (2010) seen on the 
north-east, towards DR21(OH).) Note tnat CygX-N63 is 
dominated by a single component and a weak one appears 
separated by ~0.5 km s _1 . 



3.2.2. IRAM PdBI H 13 CN spectra 

H 13 CN spectra were integrated within the 50% level of the 
peak emission and are shown here with the corresponding 
hyper fine- structure (HFS) fits for each source (Figure |3|. 
The strongest lines are observed towards CygX-N3 and 
CygX-N48, while the other cores are weaker, but all 3 hy- 
perfine components are detected. A larger line-width com- 
ponent around the main peak of the line is seen towards 
CygX-N3, -N53, but it is strongest in CygX-N63. The 
H 13 CN line-profile towards the core CygX-N63 is a single, 
broad component, which is significantly different from all 
the other fields. Two individual components can be seen 
towards CygX-N48. 



3.3. High angular-resolution integrated intensity maps 
H 13 CO + integrated intensity maps 



3.3.1. 



The integrated intensity maps obtained for the H 13 CO + 
lines are presented in Figure El a), b). For comparison con- 
tours of 3mm continuum emission from Paper I are shown. 
Note that the two maps have a similar spatial resolution of 
'. We detect strong line-emission with the PdBI in all 
cores, which resolves significant substructures of molecular 
gas compared to the IRAM 30m observations due to the 
~7 times better angular resolution. 

The integrated intensity maps of H 13 CO + including 
zero-spacings are shown in Figure [4] a). Including the zero- 
spacings allows to recover emission from larger spatial 
scales as well. The large-scale emission enhances the asym- 
metry in the emission from CygX-N48 and -N53, which 
are associated with a large-scale structure, the DR21 fila- 
ment. CygX-N53 is located towards the northern tip of the 
filament, and therefore shows an intensity distribution in 
north-south direction, which becomes stronger towards the 
main part of the filament located south in the maps. For 
field CygX-N48 the distribution of the molecular gas emis- 
sion is strongly asymmetric and drops towards the south- 
western edge of the core, close to a compact cluster of em- 
bedded young stars. On the other hand, towards the op- 
posite, north-east direction the emission gets significantly 
stronger. This shows the same morphology as the large- 
sc ale emission of H 13 CO + , which is presented in Figure 5 
Schneider et al.| fl2010| ) as well. 



m 



Figure [4] b) and c) show maps of integrated intensity 
of H 13 CO~^and H 13 CN respectively, using only interfero- 
metric data (no zero-spacings added). The distribution of 
dense molecular gas within the 0.1 pc-scale is highly non- 
uniform and shows significant substructure. The strongest 
line emission is generally observed towards the central re- 
gions of the cores, and coincides well with the 3mm con- 
tinuum peaks (indicated by grey contours). Only CygX- 
N63 is standing out because there is a prominent "hole" 
in H 13 CO + emission towards the continuum emission. The 
two lobes surrounding the continuum peak suggest that 
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H 13 CO + is highly depleted (see Section 4.1.2) and emission 



is completely absent towards the central densest regions. 

3.3.2. H 13 CN integrated intensity maps 

The emission of H 13 CN (Figure [4 c) always peaks at the 
brightest continuum sources and thus the densest parts of 
the MDCs. Comparing the distribution of the emission of 
H 13 CO + and H 13 CN, they both trace the densest parts of 
the MDCs, but are locally different. In CygX-N3 all the 
important continuum peaks are embedded in emission of 
dense molecular gas, exce pt for CygX-N3 MM4 (following 
the notations of |Paper l| ), which is outside of the H 13 CN 
emission. Also the north-western lobe strongly appearing 
in H 13 CO + , is weak in H 13 CN which may indicate chem- 
ical differences within the cores. CygX-N12 has a central 
peak embedded in an elongated, filamentary-like emission 
in continuum, which is similarly traced by weak H 13 CO + 
emission. Nevertheless in H 13 CN only the main peak of con- 
tinuum is traced. The 3mm continuum and molecular line 
emission in CygX-N48 is more diffuse and follows an east- 
west elongation. In CygX-N53 the central continuum peak 
is well traced in both molecules, but while the southern 
lobe contains two peaks in continuum and is well-traced by 
H 13 CN, it is absent in H 13 CO + . A compact 3mm continuum 
peak (probably a more evolved protostar) is found south- 
west from the central peak, which appears in H 13 CO + , but 
is not seen in H 13 CN. As mentioned before, the main peak 
of CygX-N63 is absent in H 13 CO + , but is well traced by 
strong, compact emission of H 13 CN. 

Summarizing, H 13 CN is generally more centrally peaked 
than H 13 CO + . Such differences in the distribution of the 
two lines may be explained by physically different environ- 
ments and/or abundance anomalies close to massive proto- 
stars. 



4. Analysis 

4.1. Analysis and modeling of the single-dish spectra of 
HCO+ and H 13 CO+ 

4.1.1. Line-fitting of single-dish H 13 CO + spectra 

The single-dish maps of H 13 CO + (J=l-0) show a single line 
detected over a significant extent in each core, thus we fit 
them to first order with a gaussian. We use the result of the 
fit to derive the parameters of the spectra: line intensity, the 
local rest velocity (vi sr ) and the line- width (and dispersion 
g = F ] VHM ) in each core and assume that H 13 CO + is op- 

V 81n ( 2 ) 

tically thin. The values for each core at the position of the 
PdBI phase center are listed in Table |2j 

The emission associated with cores CygX-N3 and - 
N12 have a v^ sr of +15 km s _1 , while cores CygX- 
N48 and — N53 have a v/ sr of —3.5, —4.2 km s _1 , respec- 
tively. These two cores are located in the massive filament 
associated with DR21 and their rest velocity is consistent 
with that of the filament, confirming that they are phys- 
ically associated. CygX-N63 is located southwards of the 
filament, but is associated with the south tail of the large- 
scale complex of DR2 1 and DR23 in Cygnus -X North with 
of -4.3 km s _1 ( |Schneider et al .| 
We derive velocity-dispersions 
1 km s" 



a vi s 



2006). 



Table 2. Result of line-fitting to single-dish H 13 CO + lines 
(28" angular resolution): integrated intensity, local rest ve- 
locity and line-dispersion. We compare the line-dispersions 
obtained with sin gle- dish N 2H + pointed observations to- 
wards these cores ( |Paper l[ ). 



Source 


J T mb dv 
[K] 


Vlsr 

[km / s] 


[km / s] 


(J N 2 H + 

[km / s] 


CygX-N3 


4.77 ± 0.1 


15.07 ± 0.03 


1.01 ± 0.02 


0.98 


CygX-N12 


3.30 ± 0.11 


15.27 ± 0.05 


1.23 ± 0.05 


0.89 


CygX-N48 


5.02 ± 0.14 


-3.49 ± 0.04 


1.26 ± 0.04 


1.28 


CygX-N53 


5.54 ± 0.13 


-4.18 ± 0.02 


0.9 ± 0.03 


0.76 


CygX-N63 


2.54 ± 0.09 


-4.34 ± 0.03 


0.69 ± 0.03 


0.72 



Table 3. Result of HFS line-fitting of the PdBI H 13 CN 
line, which is taken from an average within the contour of 
50% of the peak emission: integrated intensity, local rest 
velocity, line-dispersion and opacity. The same excitation 
temperature is assumed for all components. 



Source 


.[ 


/dv 


Vlsr 




&H 13 CN 


T 




[Jy/beam] 


[km/ 


s] 


[km / s] 




CygX-N3 


0.113 


± 


0.002 


14.96 ± 


0.12 


0.69 ± 0.03 


0.1 


CygX-N12 


0.076 


± 


0.002 


16.16 ± 


0.29 


0.98 ± 0.07 


0.1 


CygX-N48-l 


0.05 


± 


0.01 


-2.39 ± 


0.04 


0.53 ± 0.1 


0.15 


CygX-N48-2 


0.127 


± 


0.002 


-4.77 ± 


0.02 


0.86 ± 0.05 


0.1 


CygX-N53 


0.063 


± 


0.001 


-4.87 ± 


0.03 


1.1 ± 0.06 


0.1 


CygX-N63 


0.063 


± 


0.006 


-4.64 ± 


0.1 


2.05 ± 0.22 


0.271 



-0.1-0.2 km s -1 ) (Table [2). We note that the largest line- 
dispersions are observed towa rds the most massive core, 
CygX-N48 (see also in Section [5^2). 



4.1.2. Radiative transfer modeling 

Radiative transfer modeling of the single dish spectra of 
H 13 CO + and HCO + was performed simultaneously using 
the ID spherical, non- LTE radiative transfer code Simline 
flOssenkopfet al |200l| ). The spectra of H 13 CO+ and HCO+ 
were extracted for each MDC at the center position of the 
single-dish maps and then averaged in concentric annuli 
in steps of 14" with respect to this position. The physical 
model used for all M DCs is desc ribed in detail in Appendix 
[A] (See also Figures [Q| to | A. 7[ ). 



We used continuum data t o derive the m ass, thus the 

[2008] ) 
he den: 



the orders of 

_1 for the co res, which is similar to N2H + obser- 
vations from Paper I , though 3 of the latter are smaller (by 



density profile, and NH3 data (Wienen[2008j) to constrain 
the temperature profile and varied only the density at which 
depletion sets in and the depletion ratio for each core within 
a range of 10 5 to 10 6 cm -3 and 0-1000, respectively. The 
best fit results correspond to optically thick emission of 
HCO+ (t~8.3 to 15.4) and optically thin H 13 CO+ (r-0.3 
to 0.5) emission in all MDCs. We find that introducing de- 
pletion at a density of ~10 6 cm -3 with a depletion factor of 
50-100 gives generally good line intensities. Only within the 
coldest and densest MDC, which is CygX-N63, depletion 
had to be introduced at a lower density of ~2-10 5 cm -3 . 
With this modeling we show that for the optically thin cases 
the emission is dominated by the outer layers and not the 
regions close to the center, where depletion is high. This 
indicates that the emission is coming more from the bulk 
of material with densities between 10 4 - 10 6 cm -3 , corre- 
sponding to the inter-protostar medium, the mass reservoir 
of the core. 
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Fig. 5. a) The maps show the velocity- field of the single-dish maps of H 13 CO + , which is obtained by a single gaussian 
fit. Dashed circles show the primary beam size of the PdBL Dashed lines show the axis of the velocity gradient, while 
the arrows indicate the position angle of the gradient. The error of the gradient fitting is shown in the lower corners 
of the images, b) The maps show the peak velocit y of the PdBI+30m H 13 CO + spectra, c) Position- velocity cuts 
along the velocity gradient derived in Section 



4.2.1 



and indicated on a) from the single-dish H 13 CO + spectra-cubes. 
Contour levels go from 20% of the peak intensity by 20% levels. Dashed lines indicate the systemic velocity, d) IRAM 
30m+PdBI position velocity cuts along the same direction. Contour levels indicate 20% of the peak emission as well 
and increase by 20% for CygX-N3,-N53,-N63, while contours start at 50% of the peak for CygX-N12, -N48. 



4.2. Velocity field 

4.2.1. Global velocity field: rotation like patterns? 

The systematic shifts seen in the line-positions were stud- 
ied to understand the global kinematics of MDCs. We de- 
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Table 4. Result of least square fitting to the velocity gra- 
dient of the IRAM 30m H 13 CO + maps: velocity gradi- 
ents and positions angles. We calculated the f3 parame- 
ter (E rot /E pot ) and the specific angular momentum respec- 
tively. Note that towards CygX-N12 and -N48 no gradient 
was found. 



Source 


[km s _1 pc _1 ] 


PA 

[°] 




P 


J/M 
[km s _1 p( 




CygX-N3 


1.97 ± 0.7 


96 


30 


x 10- 4 


7.88 x 10" 


-3 


CygX-N12 


0.42 ± 0.43 


-89 


1.5 


x 10- 4 


1.68 x 10" 


-3 


CygX-N48 


0.68 ± 0.7 


46 


4.2 


x 10- 4 


5.33 x 10" 


-3 


CygX-N53 


4.17 ± 0.7 


67 


250 


x 10- 4 


20.9 x 10" 


-3 


CygX-N63 


1.17 ± 0.34 


-93 


4 


x 10- 4 


1.86 x 10" 


-3 



rive the velocity field with the gaussian line-fitting proce- 
dure to the H 13 CO + maps of the IRAM 30m telescope (see 
Figure [5]). CygX-N3, -N53 and -N63 show a homogenous, 
but clearly not isotropic velocity field with a systematic, 
axi-symmetric shift in the line position (velocity) . In CygX- 
N53 this velocity p attern follows the globa l velocity field of 
the DR21 filament flSchneider et al .||2010| ). Towards CygX- 
N12 and CygX-N48 no clear global pattern is seen, the ve- 
locity field seems to be rather random. Figure [2] b) indicates 
that in CygX-N12 a single gaussian fit may not reproduce 
well the observed line profiles due to blending of more com- 
ponents. Similarly, for CygX-N48 a second component is 
seen as a weak shoulder on the spectra, therefore a single 
gaussian fit for these two MDCs may not be representative 
of the velocity field. 

Altogether, for 3 out of 5 MDCs a rotation like veloc- 
ity pattern is visible. To derive the axis and position angle 
of the velocity shif t , we use the least square fit approach 



of 



Goodman et al. 



, ,19931 ). We applied the routine VFB^ 
( Goodman et al.|1993| ) in order to derive velocity gradients 
from the velocity field. The results are shown in Tableland 
the position angle of the resulting gradients are indicated 
in Figure [5] The derived angle corresponds well with the 
observed axi-symmetry of the velocity field. The velocity 
gradients have the same order of magnit ude of what has 



been observed for low-m ass nearby cores (Goodman et al. 



1993 



.2 km s" 1 



pc 



Caselli et al.||2002[ ), between 1.2 - 
The highest gradient is seen in CygX-N53, while in CygX 
N12 and -N48 no gradients can be derived. Among the three 
cores, where the velocity field follows a clear symmetry, the 
smallest gradient is associated with the least massive, but 
most compact core (containing only one massive protostar), 
while the largest is in CygX-N53, which is embedded in the 
DR21 filament. That itself has a very complex kinematics 
with large gradients. The velocity gradie nt of CygX-N53 
i s con sistent with previous estimates of Schneid er et al. 
(2010) indicating that this core is embedded and kinemat- 



ically linked to the filament. 



4.2.2. Small-scale velocity field 

We compare the velocity field obtained with single-dish, low 
angular-resolution data with the high angular-resolution 
observations obtained with PdBI (and zero-spacings). Since 
the spectra show several lines dispersed within a few 



3 IDL routine for performing least square fit to the velocity 
field. 



km s _1 , a single gaussian fit can not be representative of 
the velocity field. Therefore we extract and map the peak 
velocity of the spectra, which results in a velocity field trac- 
ing the bulk of motions (see Figure [5] right panel). 

At high angular-resolution the axisymmetric global ve- 
locity pattern is resolved into more substructures in ve- 
locity, indicating a much richer dynamics. Despite of the 
small-scale velocity fluctuations, in 3 MDCs the velocities 
are dominated by a large-scale gradient. CygX-N3, -N53, 
-N63 show a remarkably symmetrical global velocity field 
with systematic shifts in velocity, consistent with the single- 
dish maps. The velocity field in CygX-N12 and -N48 lack 
clear gradients and are dominated by fluctuations on small- 
scales, which make it very complex. Within CygX-N12 a 
strong outflow is observed in 12 CO (J=2-l) (Bontemps et 
al., in prep), which may also be apparent in the velocity 
field map as indicated in Figure [5] Towards CygX-N48 the 
velocity field is complex without any clear trends of system- 
atic motions. The complexity of such velocity field towards 
CygX-N12 and CygX-N48 is due to the presence of several 
spectral components, which are also seen on Figure [2] 

4.3. Position-velocity diagrams 

Position-velocity cuts were taken along a cut through the 
continuum peaks and perpendicular to the axis derived 
from the fit to the velocity fields. In Figure [5] c) contours of 
the position-velocity maps of the single-dish observations 
are shown. The velocity gradients towards CygX-N3, -N53 
and -N63 are clearly seen. For CygX-N12 and -N48 a large 
dispersion in velocity becomes obvious, which explains the 
complexity of the velocity field maps. 

Figure [5] d) shows the position- velocity maps along the 
same cuts from the high-resolution IRAM 30m+PdBI data. 
In all cases we resolve much richer structures confirming 
the existence of several line components and gradients. For 
CygX-N3, -N53 and -N63 these velocity cuts illustrate per- 
haps best the complexity of velocity features obtained with 
the high angular-resolution data. In CygX-N3 the pv-map 
shows a clumpy structure with several peaks showing a gra- 
dient from ~14.5-16 km s _1 . The peaks correspond to in- 
dividual structures at different positions, but we point out 
that they must be physically linked, being embedded in the 
same flow of gas and consistent with the velocity gradient 
seen in the single-dish data. In CygX-N53 and -N63 the ve- 
locity shift is smooth (less clumpy) and goes from ~-5.5 to 
-3.5 km s _1 and from -3.5 to -5 km s _1 , respectively. 

For CygX-N12 and -N48 we resolve for the first time in- 
dications for velocity gradients. On these two plots (Figure 
[5] d, middle panels) we increased the contour levels to show 
only the 50% level of the peak emission, since there is 
emission coming from broad line- wings due to outflows. 
For CygX-N12 the gradient goes from 17 km s _1 down to 
14 km s _1 , while CygX-N48 shows two components: one 
dispersed between -2 and -4 km s _1 , the other component 
at ~-5.5 km s _1 . 

Such velocity gradients, as seen for CygX-N3, -N53 and 
-N63, can be interpreted in the frame of rotating structures. 
In case of solid body rotation, velocity scales linearly with 
the position, which could be consistent with the pv-maps for 
CygX-N3, CygX-N53 and CygX-N63. For the other fields 
solid body rotation is unlikely due to the complexity of the 
individual velocity components. The origin of the velocity 
pattern is discussed in Section [5.2.2 
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4.4. H^CN spectra 

We performed HFS line-fitting to the extracted H 13 CN 
spectra for each core and summarize the results in Table 
[3j For CygX-N3, -N12, -N53 and -N63 we obtain a good fit 
with only a single line, but we indeed identify two compo- 
nents for CygX-N48 at -4.77 km s" 1 and at -2.39 km s" 1 . 
The line velocities of H 13 CN are in good agreement with the 
velocity of the H 13 CO+ line (±0.2 km s -1 ) towards CygX- 
N3, and -N63, while for CygX-N12 and -N53 they are off 
by ~0.8 km s _1 compared to the H 13 CO + line. The rea- 
son behind this is not clear, but may indicate that there are 
chemical differences on the projected size-scales of ~ 10.000 
AU and these two tracers may not trace the same gas. The 
HFS fitting in Figure [3] suggests that the assumption of 
LTE excitation condition is satisfactory for all MDCs. 

In Section [3.2.2| we noted already the "bump" in the 
spectra, peaking at the position of the main line with a large 
plateau on the two sides. We further investigated the possi- 
ble origin of this feature by searching for complex molecules 
at these frequencies but found no lines that could result in 
such a strong contribution. This is the most prominent fea- 
ture in CygX-N63, where the line-dispersion is surprisingly 
large (~2 km s _1 ) compared to the other MDC s. We discuss 
in detail the possible origin of that in Section [5.1.2 



5. Discussion 

From a theoretical point of view, there is either a high level 
of turbulence in the MDCs, capable to provide a significant 
support against gravity to allow a quasi-static evolution 
towards the monolithic collapse of the cores, or the cores 
are not in equilibrium and their formation and evolution is 
mostly driven by dynamical processes. Specific and different 
kinematics is expected in the two scenarios (e.g. supersonic 
micro-turbulence versus complex, but organized flows). 

W e firs t discuss the origin of the observed line emission 
(Sect. |5.1| ), then we analyze the physical origin of the ve- 
locity dis pers ion at large (~0.1 pc) and small scales (~0.03 
pc, Sect. 5.2) to ultimately discuss the amount of micro- 



turbulence in the cores which could stabilize and regulate 
their evolution. We finally review the obtained resu lts i n 
the context of the two scenarios in competition (Sect. [573]) . 



5.1. Molecular lines to trace high-density material 

5.1.1. H 13 CO + traces the inter- protostar medium 

H 13 CO + has often been detected around low-mass pro- 
tostars following an elongated structure perpendicular to 
their outflow. Having a high critical density, it is acknowl- 
edge d as a good tracer of low-mass pro tostellar envelopes 
(e.g. |Saito et aLl|2001| |j0rgensen et al.||2004| . The average 
density in o ur sample of MDCs is high (>10 cm -3 ) ( |Motte| 
| et al. | 2007 ), while the peak density in protostars estimated 
" is reaching up to -10 7 -10 8 cm . These MDCs 



m 



Paper I 



lack significant emission at infrared wa velengths, and re- 
cent NH3 measurements ( Wienen||2008 ) confirm that they 
are cold with kinetic temperatures between 18-28 K. Such 
high densities at low temperatures should lead to a high 
level of depletion of virtually all C-bearing molecules as a 
consequence of CO depletion onto grains (e.g. Tafalla et al. 
2002} |Bergin fc Tafalla|2007j ). On the other hand, the des- 
orption temperature for CO is of the order of 20 K, i.e. 



similar to the observed temperatures here, considering also 
that if protostars have formed already, they might heat-up 
their cold environment, releasing and restarting the forma- 
tion of C-containing molecules. The precise level of deple- 
tion is therefore uncertain and the abundances of different 
molecules (here H 13 CO + and H 13 CN) as a function of den- 
sity and temperature is even more difficult to predict. 
From the radiative tran sfer m odeli ng (Section |4.1.2 



A. 7) we conclude that 



Appendix [A] and in Figures A. 3 
H 13 CO + is always optically thin and the emission is dom- 
inated by the gas between the protostars and can not 
probe their highest densities due to significant depletion 
in their central regions. We find indications of the highest 
depletion towards CygX-N63, where in its central regions 
a significantly higher depletion ratio and lower depletion 
density was used to reproduce the correct line intensities. 
Summarizing this picture, we base the following discussion 
on considering that the emission of H 13 CO + comes from the 
dense material within the MDC, but rather representing 
the bulk of material in which the protostars are embedded, 
than their inner envelope as for low-mass protostars. 



5.1.2. H 13 CN, a better tracer of the highest density regions? 



Modeling of e.g. Lee et al. (2004) has shown that HCN may 
also be affected by depletion, but to a somewhat lowe r level 
than HCO + . In surveys for low-mass pre-stellar cores [Sohnl 
|et al.| (2007) observed that HCN traces better the inner part 
of the cores than HCO + and shows a similar distribution 
as N2H + , which is known to trace the inner (dense, cold) 
part of the cores ( Tafalla et al. 20 06; Di Francesco et al. 
2007). Nevertheless, |Lee et al. (2004) also pointed out that 
HCN is part of a complex chemical network, whereas HCO + 
depends directly on the CO abundance. Based on the mor- 
phology of the distribution of emission presented in Figure 
[IJ we show here that H 13 CN has a better correspondence 
with the central higher density regions (continuum peaks), 
than H 13 CO + . The most prominent example is CygX-N63, 
where H 13 CO + is completely depleted while H 13 CN peaks 
at the strong continuum source. The origin of this may be 
due to two reasons: H 13 CN may be simply less depleted - 
as it has already been pointed out in examples of low-mass 
cores, and/or H 13 CN may be affected by different chem- 
istry, with for instance overabundance due to strong shocks 
emerging from the close vicinity of the protostars. 



5.2. Micro-turbulence and bulk motions within the MDCs 

Analysis of line-widths reveals the contribution of ther- 
mal versus non thermal motions. The non thermal ve- 
locity dispersion is usually interpreted as isotropic turbu- 
lence down to small scales where it is often referred to as 
micro-turbulence. However it may have contribution from 
bulk motions, such as rotation, infall and coherent flows or 
streamings. It is of high importance to measure the true 
level of micro-turbulence since this leads to a local support 
of similar nature as thermal motions and thus contribute 
to the global support of the cores against gravity. 

We first discuss the global vel ocity d ispersion as com- 
pared to previous studies (Sect. |5.2.1 ). The analysis of 
the single dish observations showed that the line widths 
of H 13 CO + for most of the observed MDCs have con- 
tributions of systematic, global motions which are not 
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micro -turbu lence, namely rotation (Sect. 5.2.2) and infall 
(Sect. 5.2.3). The high resolution PdBI data confirm these 
motions but, in addition, disentan gle sev eral velocity com- 
ponents at small scales (see Sect. 5.2.4 below). Finally, it 
is only after these different contributions to the line- widths 
are recognized, that we will be able to discuss the possible 
remai ning level of micro-turbulence to support the cores 
(Sect. K3). 



5.2.1. Global velocity dispersion and virial equilibrium 



2.0 
1 .5 
1 .0 
0.5 
0.0 
-0.5 
-1.0 



A one-dimensional thermal velocity 
H 13 CO+and H 13 CN gives a T = 



for 



dispersion 
where ks 
li n is the 



-1 .5 




3ive cores (Caselli et al.) * 
H 13 C0 + Cygnus-X MDCs • 
N 2 H + Cygnus-X MDCs + 
RDCs (Plume et al) X 



IS 



-2.0 



-1 .5 



UN rn H 

the Boltzmann constant, T is the temperature 
molecular weiec ht equal to 30 and 28 for H 13 CO+and 
H 13 CN respecti vely and trh is th e mass of a hydrogen 
atom. Following |Motte et a l. (2007) and Paper I we assume 

a temperature of 20 K, which gives a thermal velocity- ture ( |Caselli fc Myers 
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Log R [pc] 



0.5 



Fig. 6. We plot line-widths (FWHM) of the Cygnus-X 
MDCs (red large dots and blue crosses) versus the core size. 
Black symbols indicate line-widths taken from the litera- 



dispersion of 0.074 and 0.077 km s" 1 for H 13 CO + and 
H 13 CN respectively, while the observed dispersion is 
~1 km s _1 (Tabled). This is much (~15x) larger than the 
thermal velocity-dispersion, indicating that non-thermal 
motions significantly contribute to the line broadening. 

Figure [6] shows the line- widths of the single dish 
H 13 CO + (this work) and N 2 H + lines (from |Paper l[ ) for 
the 5 MDCs in comparison with sa mples of single-di sh ob- 
servations of massive cores in Orion (Ca selli fc M yers 1995) 
and Infrared Dark Clouds (IRDCs) ^Plume et al ||1997p . 
The Cygnus X M DCs appear to be interm ediate between 



the Orion cores of [Caselli fc Myers] (|1995 ) and the larger 
P 



distance IRDCs of Plume et al. 



We use Eq. 18 and 20 wit h t 



(1997JT 

;he fiducial model assump- 
tions of McKee & Tan ( 2003 ) to predict the required tur- 
bulent velocity dispersions at a corresponding core radius 
to keep the MDCs in equilibrium: 



1.27 



V30M^y 



1/4 



0.057 



\30MrJ 



a/4 



1/2 



-1/2 



(1) 



(2) 



where the parameter m*f corresponds to the final mass 
of the star and E c j denotes the average clump surface den- 
sity. We calculate the turbule nt dispersion and the core size 
adopting £ cZ = lg cm" 2 (like |McKee fc Tan|2003| > for a 10 
to 60 M final star. This adopted surface density is com- 
parable or at least a fact or of 2 l ower than the observed 
surface density of MDCs (Paper I). 

Our measurements at the scale of the MDCs (0.1 pc 
scale), which themselves indicate upper limits for the turbu- 
lent support (due to other contributions to the line- width), 
are s maller than the model predictions of McKee & Tan 
(2003). In fact the maximum level of micro turbulence, 
which is observed, would only support cores with surface 
density of a factor of 3 lower than E c j = lg cm -2 (see 
dashed line in Figure |6|. 

5.2.2. Angular momentum and rotation of the MDCs 

The velocity field derived from H 13 CO + maps at the size- 
scale of the cores (~0.1 pc) indicates that non-thermal, 
organized motions such as rotation also contribute to the 



like line- widths size-re' 



1995). Black line indicates a Larson 
ation^while the dashed-ditted line is 



adopted from |Plume et al.| ( |1997[ ) representing their best ro- 
bust fit to their observed line- width size-relation of IRDCs. 
Grey line indicates the fiducial model of McKee & Tan 
( 2003 ) with a canonical surface density of ~ lg cm -2 and a 
final stellar mass of 10 to 60 M . Dashed grey line indicates 
the same model with a clump surface density 3 times lower. 
A small box indicates the zoom presented in Figure 10 



observ ed line- widths. According to Burkert & Bodenheimer 
(2000) if the cores were dominated by purely isotropic 
micro-turbulence, the velocity field should be completely 
random. From Figure [5] it is clear that CygX-N3, -N53 
and -N63 show indications for organized motions, while for 
CygX-N12 and N48 this is less obvious. For the before men- 
tioned MDCs the velocity maps are axisymmetric and the 
axis of the symmetry is close to be aligned with the mini- 
filament containing the protostellar fragments. As a first 
order assumption we interpret this phenomenon as rota- 
tion, which is in fact a natural consequence of preserving 
the initial angular momentum of the col lapsing cloud. We 



note that recent numerical simulations by Dib et al. (2010a 
do report a variety of dynamical patterns for cores formed 
in a turbulent, magnetized, and self- gravitating molecular 
clouds ranging from easily recognizable rotational features 
such as the ones observed in CygX-N3,-N53, and -N63 to 
more complex features, such as the ones seen in CygX-N12 
and -N48. 

To further discuss the observed angular momentum we 
rely on the systematically derived velocity gradients. Our 
values of 1.2-4.2 km s _1 pc _1 are consistent with those re- 
ported in the literature of both low and high-mass objects. 
We also calculated the (3 parameter (see Table |4|, which 
shows the ratio between the rotational and potential en- 
ergy (E rot /E pot ). Values between 2 x 10 ~ 3 and 1.4 were 
report ed in a sample of low-mass cores by |Goodman et al.| 
fll993| ). Our f3 values are low (4 x 10~ 4 - 2.5 x 10~ z ) (see 
Figure [7] b) indicating that gravity significantly dominates 
over the rotational energy, especially in CygX-N3 and - 
N63. The most extreme example is the most compact core, 
CygX-N63 (/? ~ 4xl0 -4 ), which hosts a single or a close 
binary massive protostar. A priori one could expect this 
core to behave similarly to low-mass cores, where rotat- 
ing envelopes around Class objects are reported. Its low 
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Fig. 7. a) Plot of size versus velocity gradient comp aring 
a sample of low-mass cores (Goodman et al. 1993) with 
the Cygnus-X MDC sample. Red symbols represent the 
Cygnus-X MDCs, whi le black stars show the sample of 
oodman et al. (11993) . Black line shows the best fit from 
Goodman et al. (1993) with <S = io-°- 3 ± - 2 R-°- 4 ± - 2 . 
b) Plot of p versus size relation shown for the same sample, c) 
Plot of the specific angula r momentum versus siz e. Black line 
shows the best-fit line from Goodman et al. (1993) with J/M = 



10 



-0.7±0.2p£l.6±0.2 



z 0. 



velocity gradient corresponds to a very low specific angu- 
lar momentum (F igure [7| c) fitting on th e trend found for 
low mass cores by |Goodman et al.| ( |1993[ ) . Nevertheless the 
other MDCs have higher specific angular momentum. We 
suggest that the origin of this diversity is related to the 
small-scale dynamics, and that the origin of the velocity 
field may be differen t for CygX-N63 and the other cores 
(see in Section 5.2.4). 

In contrast, towards CygX-N12 and -N48 no clear 
rotation-like gradient are evident in the large-scale (~0.1 
pc) velocity field maps with complex velocity patterns (see 
left panels of Figure [5]). This suggests that these 2 MDCs 
do not have a large resulting rotation in projection on the 
sky. Their main rotation axis may be close to the line of 
sight or these MDCs have a truly smaller angular momen- 
tum than the others. In the high angular-resolution maps, 
indeed some velocity gradients are observed at small scales 
with similar velocity drifts like in the other MDCs, indicat- 
ing that the angular momentum may be redistributed at 
smaller scales. 



The position-velocity cuts along the derived velocity 
gradient (Figure [5] right panel) show a centrally peaked 
elongated pattern for CygX-N3, -N53, -N63, which may cor- 
respond to large rotational toroids up to 40" (68.000 AU), 
45" (76.500 AU) and 35" (60.000 AU) respectively, which 
may serve as a common and coherent reservoir of mate- 
rial for forming protostars. Such rotating structures with 
similar size-scales were observed towards low-mass Class 



0/Cla ss I objects (e.g. Goodman et al. 1993 Caselli et al. 



2002 ) and towards massive star- forming cores as well (e.g . 



Beltran et aL]|2004j |Cesaroni|[2005l |Fallscheer et al.||2009 ) 

but the nature and origin of these rotating toroids are not 
well established yet, especially in cores hosting the earliest 
stages of massive star-formation. 



5.2.3. Importance of global infall 

The Cygnus-X MDCs have larger masses for the same sizes 
and slightly larger velocity dispersions compared to the 



Orion cores. The virial-masses are determined in Pa per I 
and comparing the virial-mass to mass ratios (1.16, 0.87, 
1.06, 0.75, 0.51 for CygX-N3,-N12,-N48,-N53,-N63 respec- 
tively), we find that most of the Cygnus-X MDCs are grav- 
itationally unstable and are therefore expected to be in 
global collapse. Four out of the five MDCs indeed display a 
typical blue-shifted asymmetric HCO + line pr ofiles suggest- 
ing infall motions. As showed in Sect. |4.1.2[ the line pro- 
files of both HCO + and H 13 CO + can be reproduced with 
radiative transfer modeling with typical infall velocities be- 
tween 0.1-0.6 km s _1 . This indicates that in addition to the 
rotation-like organized motions discussed above, the global 
velocity dispersions in the MDCs have also a contribution 
from infall motions. The strongest infall is seen towards 
the most massive core CygX-N48 (~0.6 km s _1 ) and is 
similarly large towards CygX-N53, i.e. the 2 MDCs l ocated 



in the DR21 filament studied by |Schneider et ah] ( [2010 ) . 
CygX-N3 shows a red-shifted asymmetric profile, which is 
indicative of expanding motions. However the stronger peak 
observed at ~16.5 km s _1 is coinciding with a second ve- 
locity component identified only from th e hig h-resolution 
spectra map in Section p] (see also Table B.l) and there- 



fore we argue that this peak has contribution from an ad- 
ditional velocity component, rather than expansion. HCO + 
emission towards CygX-N63 is significantly weaker than in 
the other cores and radiative transfer models show that de- 
pletion is the strongest in this core. Nevertheless the HCO + 
line shows a slight, blue-shifted shoulder, which also points 
to infall-motions. Our best radiative transfer model gives 
an optically thick line for HCO + with an in fall speed of 
~0.1 km s _1 . The high depletion, as shown by|Rawlings fc| 



| Yates ( 2001| ) based on numerical simulations, might explain 
this rather weak indication of infall. 



5.2.4. Coherent flows seen at high spatial resolution 

The high spatial resolution PdBI velocity maps in Figure [5] 
(second column of panels) display velocity fields which are 
similar to what is seen in the purely single dish data (first 
column of panels in Figure [5| in particular the gl obal 
rotation-like gradients discussed above in Sect. |5.2.2| On 
the other hand these global motions are found to have a 
more complex morphology at small scales with stronger lo- 
cal gradients, and more importantly with the H 13 CO + line 
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Fig. 8. Map of H 13 CO + emission obtained with the PdBI and zero-spacings. The grey scale shows integrated intensity 
around the v^ sr , between 14.5-15.5 km s -1 . Blue contours correspond to integrated intensity between 13.5-14.5 km s _1 
and red contours show integrated intensity between 15.5-16.5 km s _1 . Contours go from 5xrms noise and increase in 
steps of 3xrms noise. The black dashed contour shows the 30% of the peak intensity of the 3mm continuum map. Fully 
sampled spectra are extracted in two cuts following two filamentary structures, where the position corresponding to each 
spectra are shown in the upper right corner. The integration ranges from the blue shifted component are in dark grey and 
white grey shows the red component. The vi sr is shown by a dashed grey line. The right panels show position-velocity 
diagrams for both cuts respectively. 



profiles, which are not single gaussians and clearly split into 
several individual velocity components. 

These individual components may correspond to coher- 
ent features in velocity. To trace the bulk emission from 
the mass reservoir we integrate around the rest velocity of 
each MDC and then separate the emission from the other 
components by integrating in the blue and the red side of 
the velocity range. This representation allows to study the 
distribution of gas associated with the bulk emission and 
traces the morphology of the potential coherent velocity 
features of these individual components. 

Figure [8] shows CygX-N3, where the bulk emission at a 
rest velocity of 15 km s _1 (grey scale) is coinciding with the 
main continuum peaks (dashed black contours) and shows 
an elongated, filamentary structure centered on it as well 
as hints for sub-filaments in the north-east and south-west 
edge, perpendicular to it. The blue and red-shifted spectral 
components show a bipolar structure indicative of veloc- 
ity gradients or shears along these sub- filaments. In order 
to further investigate the velocity structure of these sub- 



filaments, we show posit ion- velocity cuts along them. The 
upper cut shows a global smooth shift in velocity, but with 
hints for two sharp transitions or jumps at offset 16", 3" and 
6", 10". The lower cut shows also a strong red component 
with a smooth shift in velocity and at ~15 km s _1 a sharp 
change in velocity between offsets —4", —5" and — 5",— 4". 
Such velocity shears may be similar to the vel ocity jump 
revealed in NGC2264-C by 



Peretto et al. 



(2006) 



We show spectra across the cut on the sub-filaments. 
In the upper one, until offset 16", 3" the blue component 
dominates and seems to stay at the same velocity, then the 
spectra suddenly get broader peaking at the rest velocity. 
After 8", 9" the red component appears and becomes more 
dominant at 0", 15", even if the separation between the 
main and the red component is less clear. The position of 
the red component is also compatible with a fixed velocity 
along the sub-filament. The separation between the blue 
and the bulk emission is ~0.75 km s _1 and ~0.55 km s _1 
between the red and the bulk emission. The lower cut shows 
both the blue and red spectral components from -l",-8" to 
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Fig. 9. CygX-N63 presented similarly as Figure[8j Integration ranges for the grey scale is -4.0 to -4.7 km s _1 , red contours 
are integrated between -2.7 to -4.0 km s _1 , blue contours are integrated between -4.7 to -5.7 km s _1 . Contour levels go 
from 5xrms noise by steps of 3xrms noise. 



-5 // ',4 and a velocity jump from the blue to the main compo- 
nent is seen at -6 // ,-2 ,/ offset. It is apparent, that the velocity 
separations are much larger in the lower, than in the upper 
cut up to a total velocity difference of ~2.5 km s _1 between 
the red and the blue component in the lower cut. Towards 
the continuum peak the spectral components merge into 
a single, broader line one suggesting physical interaction 
between the two components. Summarizing, velocity shifts 
along organized gas structures are seen but a more careful 
analysis of the line profiles suggests rather local shears than 
smooth velocity gradients between the individual velocity 
components. Interestingly these velocity jumps correspond 
to weak continuum sources (see Figure pT|), which may be 
at an early stage of their formation. 

In contrast CygX-N63 is shown in Figure [9j where the 
bulk emission at a rest velocity of —4.3km s - is similarly 
centered close to the continuum peaks and shows hints of 
a north-west to south-east organized elongated structure. 
The red and blue shifted emission shows a bipolar structure 
perpendicular to the main filament. The position- velocity 
cuts along these structures show though a much smoother 
shift in velocity or a weaker velocity shear. Cuts of spec- 
tra also show that the line-profiles are less complex as for 
CygX-N3 and could be compatible with a single line com- 
ponent, which shows a shift in velocity of ± 0.6 km s _1 . 



Altogether, CygX-N63 seems to be more ordered, than the 
example above. 

Following this systematic way, we find similar trends 
in the other fields, but due to a larger complexity of indi- 
vidual velocity components they are less clear for CygX- 



N12, -N48 and -N53 respectively (see Figures C.l. C.2 



C.3 in Appendix O. Cuts of spectra and position- velocity 
maps are presented therefore centered on the continuum 
peaks and perpendicular to each other in order to give a 
glimpse of the spatial distribution of the velocity compo- 
nents. Towards CygX-N12 the main component is domi- 
nated by spherical extended emission with a large num- 
ber of individual components and no clear morphology can 
be found for the blue and red components. As the spec- 
tra show, the two components are both visible over a large 
extent of the map with a velocity of ~14.3 km s _1 and 
^15.7 km s _1 . These lines show either a shift in velocity 
over ~±1 km s _1 with respect to the systemic velocity of 
15.2 km s _1 or these may be additional components. 

Similarly, towards CygX-N48 several line components 
are detected within a velocity range of 6 km s _1 correspond- 
ing to the most complex dynamics in the sample. Two main 
components can still be separated, one around -5.5 km s _1 
and one around -3.8 km s _1 , both showing shifts in ve- 
locity and additional components appear as well indicating 
complex kinematics. 
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Fig. 10. We plot line- widths of the MDCs taken from 
single-dish and high angular- resolution data. Th e corre- 
sponding sizes are taken from Motte et al. ([2007) for the 
single-dish data, while for the high angular-resolution data 
sizes show the synthesized beam of the PdBI. The Larson- 
scaling is indicated by a dotted line. G rey line indicates the 
fiducial model of McKee & Tan ( 2003[ ) with a canonical sur- 
face density of ~ lg cm~ z and a final stellar mass of 10-60 
Mq. Different colors indicate the different tracers observed 
with only single-dish or with the interferometer. 



In CygX-N53 the bulk emission at a rest velocity of 
—4.2 km s _1 shows a north-south elongation, however the 
determination of the rest velocity may be biased due to the 
surrounding large-scale filament, which may have a redder 
bulk emission. Bipolar blue and red-shifted emission is seen 
close to the continuum peak, but the tendency is less clear 
as for CygX-N3. 

We conclude that, except perhaps CygX-N63, all MDCs 
of the sample display a mixture of several velocity compo- 
nents separated by up to several km s _1 , and show drifts 
and discontinuities in velocity. The spatial distributions of 
these components suggest that they are associated with the 
MDCs and indicate that they do not show up as only indi- 
vidual fragments which could be associated with the dust 
continuum fragments discussed in Paper I. Instead, they 
appear as filamentary structures, showing systematic ve- 
locity differences leading to local shears. For three of the 
five MDCs, these velocity features averaged at the scale of 
the MDC result in a rotation-like motion, and for four of 
them, there is also an indication of global infall. This is like 
if the global motions (rotation and infall) split into small 
scale individual streams or flow s when observed at high spa- 
tial resolution (see Sect. 5.5.1 for a further discussion of the 
nature of these flows). 



velocity dispersions of the individual components decrease 
from ~1 km s _1 at ~0.1 pc-scale to ~0.6 km s _1 at ~0.03 
pc-scale. 

Velocity dispersions towards the central regions (3mm 
continuum peaks, indicated by contours in Figures [8|9| 
were extracted as well. Interestingly, the average velocity 
dispersion is larger ~ 0.86 km s _1 towards the center. 

These velocity disper sions are close to the predictions of 
the McKee & Tan (2003) model. However the line-widths in 
the central regions might not be representative of the bulk 
of the central dense material for two reasons. For H 13 CO + , 
the high level of depletion makes that only the outskirts are 
traced, which could have a larger velocity dispersion and 
H 13 CN may trace only higher temperature gas influenced 
by stellar feedbacks (outflow shocks). Alternatively, local 
turbulent s upport could be compatible with the McK ee fc 
Tan ( |2003| ) model at these small scales in the central re- 
gions. 

On the other hand, we point out that the turbulent line- 
dispersions of the individual velocity components, which are 
generally found off-center, represent more the pre-collapse 
gas and are smaller than the model values. They are found 
to fall on the Larson-scaling. 



5.4. Are massive-star forming cores in equilibrium? 

The fundamental question is if MDCs are in a quasi-static 
equilibrium or are go verned by dynamic pr ocesses. In the 
scenario proposed by |McKee fc Tan ([2003 ) massive cores 
should evolve in a quasi-static equilibrium, where a high 
effective turbulent pressure balances gravity. Once grav- 
ity becomes stronger than the effective pressure within the 
core, it collapses monolit Ideally to form a single or a bi- 
nary object. During the quasi-static evolutionary stages a 
homogenous medium should develop a central concentra- 
tion of matter while it is stabilized by its internal micro- 
turbulence. Comparing our single-dish observations of op- 
tically thin emission (H 13 CO + ), we confirm a centrally 
peaked distribution of emission in a homogenous medium. 

In the sections above we systematically searched for 
such a support, that can be micro-turbulence, magnetic 
fields and rotation. Neglecting magnetic fields as it is not 
known, we found that: a) line-dispersions on a 0.1 pc-scale 
tend to show that micro-turbulence alone is not sufficient 
to stabilize MDCs using the turbulent-core model; b) line- 
dispersions at small-scales indicate that the initial condi- 
tions may be even less turbulent; c) rotation like ordered 
bulk-motions are not sufficient either to stabilize MDCs 



(Section 5.2.2). Altogether this gives a coherent view of 
MDCs being out of equilibrium state, as indicated also by 
global infall motions and supersonic flows. 



5.3. Importance of the turbulent support at small-scales 

A detailed understanding of global organized motions both 
at low and high angular-resolution is required to recognize 
the r ole of support by micro-turbulence at small-scales. In 
Sect. |5.2.1| it was pointed out that turbulence within the 
MDCs at global scales (~0.1 pc) tends to be smaller than 
required by the turbulent core model. Here we can further 
investigate this trend at smaller scales, where the level of 
micro-turbulence is traced by the line-dispersion of the in- 
dividual spectral components recognized in the previous 
section. Figure [To] shows these results, where the average 



5.5. Dynamical processes in the MDCs 
5.5.1. Evidences for converging flows 

Several indications are found for dynamical processes act- 
ing at small-scales. The individual velocity components 
seen in the high resolution H 13 C O + lin es are best under- 
stood as individual flows (Sect. |5.2.4 ). Since the MDCs 
are in probable global collapse, these flows are expected 
to stream inward. They should then ultimately interact in 
the central regions to build up high density cores at the 
stagnation points, creating the seeds for the protostellar 
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objects. We estimated the crossing times of these conver- 
gent flows by dividing the characteristic size of the core 
(on average R core =0.1 pc) with the characteristic veloc- 
ity shift measured by the largest separation of the indi 



vidual spectral components (<5v) discussed in Section 5.2.4 
10 5 R < 



9.£ 



[pc] [km s 1 ] ' 



In Table 5 we list these pa- 



rameters and indicate the free-fall times \jff) from Paper I 
(For CygX-N63 the crossing time estimate is based on the 
velocity dispersion of the spectra, since the individual ve- 
locity components are less clear.) To calculate the crossing 
times we adopted a correction for a possible inclination ef- 
fect assuming an average angle to the line-of-sight of ~57.3° 
(which is the average angle assuming random distribution 
of orientation angles). We give also the range of r cross for 
variations of the angle by ~20° to give a hint of the effect 
of the inclination uncertainty. 

We find that all crossing times are comparable to the 
free-fall times, which indicates that these flows are the main 
driver for the evolution of MDCs, and can clearly play the 
major role in building up material of which massive proto- 
stars form. 

The origin of these small scale flows can just be a natural 
consequence of the large scale turb ulent nature of the par- 
ent clumps, and clouds. At pc-scale [Schneider et al. (2010) 
showed that massive structures, like the DR21 filament it- 
self, can be formed by dynamical processes and converging 
flows, pointing out that at sub-parsec scales dynamical pro- 
cesses are important as well, leading to a complex view of a 
hierarchical dynamics from pc-scale down to the protostar 
scale. 



5.5.2. Importance of competitive accretion 



I n a broad view of t he competitive accretion scenario 
Bonnell & Bate 2006), the above discussed high level of 
ynamics necessary leads to competition for mass between 
the protostars. Even if such small-scale flows converging to 
the central parts of MDCs were revealed in the previous sec- 
tions, it is not clear whether these flows continue to build 
up the mass reservoirs of the protostars. The time-scale to 
have a new prestellar/protostellar seed formed by conver- 
gent flows is given by the typical velocity shears observed 
in the most dynamical regions of the MDCs at a scale of 
the or der of the separation between the protostars (~5000 
AU in Paper I). The largest velocity shifts observed over 
this projected distance are of the order of 3 km s _1 which 



lead to a typical time-scale of ~ 10 yr , a factor of a few 
smaller than the crossing time at the scale of the MDCs. 
This means that during a crossing time of the MDCs, sev- 
eral protostars could form by the convergence of several 
small scale flows. This time-scale of ~ 10 4 yr is also only 
slightly larger than the typical free-fall time derived for the 
continuum fragments reported in Paper I (on average of 
4000 yr). If the protostellar cores collapse in a few free-fall 
times like in low-mass star-forming regions, the formation 
time-scale could be of the same order of magnitude than 
the collapsing time, and the mass input from flows during 
collapse can be significant. This therefore points to the pos- 
sibility of a significant competitive accretion in the MDCs 
of Cygnus X. 



Table 5. Parameters corresponding to the identified flows. 
The velocity difference indicates the relative velocity of the 
flows on the global scale of the MDCs, which was used to 
calculate the crossing time-scales. 



Source 


5v 


Tcross 


T ff 




[km s -1 ] 


xlO 4 [yr] 


xlO 4 [yr] 


CygX-N3 


2.5 


2.4 - 3.8 


6.0 


CygX-N12 


3 


2.0 - 3.2 


5.9 


CygX-N48 


3.5 


1.7 - 2.7 


6.1 


CygX-N53 


1.1 


5.4 - 8.7 


7.6 


CygX-N63 


0.5 


11.8 - 19.1 


3.4 




Aof (arcsecond) 

Fig. 11. The velocity field map for CygX-N3 is shown with 
overlays of integrated intensity (contour levels are from 20 
to 50cr by 5<j. White contours show the 3mm continuum 
emission. Arrows indicate the shears associated with weak 
continuum emission. 



5.5.3. Origin of the diversity between the 5 MDCs 



From the fragmentation properties (Paper I) and a kine- 
matical point of view (this work) our sample of MDCs is 
diverse. This diversity can reflect either different evolution- 
ary stages or different initial conditions for the MDCs. 

In CygX-N12, -N53, -N63 the gravitational potential is 
dominated by a few, but massive fragments in the center, 
corresponding to more than 30% of the total mass. For 
these cases, the self-gravity of the MDC is dominated by 
these few fragments which could then collect most of the 
remaining flowing mass if the convergent flows are still ac- 
tive. In CygX-N3 and -N48, on the other hand, the mass 
in the already formed fragments is a smaller fraction of 
the total mass, therefore the local gravitational field is not 
dominated by the observed fragments, but more by the self- 
gravity of the mass reservoir of the MDC. These MDCs are 
either intrinsically different from CygX-N12, -N53, and - 
N63 or they are younger and will evolve to similar MDCs 
dominated by a few massive fragments. The high resolu- 
tion velocity field in H 13 CO + should indicate whether or 
not the fragments are still collecting mass from their par- 
ent MDCs, and therefore may help to discriminate the two 
possibilities. 

In CygX-N3, the 3 most massive fragments are lo- 
cated on a significant shear in velocity which indicate 
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that they can still gain mass from the MDC reservoir. 
At the same time three shears associated with significant 
H 13 CO + emission peaks are also visible and could therefore 
correspond to newly forming protostars (all three corre- 
spond also to weak continuum peaks as indicated by arrows 
in Figure 11 ). The relative intensity of the shears and asso- 
ciated H 1 ^UO + peaks suggest that a large number of pro- 
tostellar seeds collect mass from the MDC mass-reservoir 
leading to the formation of a whole cluster of protostars. 

In CygX-N12, several velocity components imply a com- 
plex kinematics. Since 70 % of the total mass of the MDC is 
still available, there is a potential for CygX-N12 to further 
evolve towards a state similar to CygX-N53 with two main 
massive fragments in the central regions. 

In CygX-N48, the observed shears are intense and are 
located in the central regions close to, but not necessary 
exactly associated with the main central fragments. This is 
like if the first fragments are not yet clearly decoupled in 
the central, very active region. This is therefore a good can- 
didate MDC to be in an early phase of evolution, perhaps 
leading later on to a MDC similar as CygX-N53. 

In CygX-N53, most of the total mass is contained in 
two main, central fragments. No clear shears are observed 
in the central regions. Some indications of shears are seen 
in the southwest outskirts of the MDCs, though they are 
less convincingly associated with the MDC (this MDC is 
located in the north part of the massive DR21 filament). 
This MDC is therefore probably evolved with no significant 
flowing gas present in the central regions. 

Finally, CygX-N63 is a special case with only one cen- 
tral peak of emission in the continuum such as if it were 
not sub-fragmented, and with a global collapse towards this 
single object. The kinematics shows a rotation-like gradi- 
ent which may correspond to a large scale shear, revealing 
that the object still collects mass. Interestingly, besides this 
smooth shift or a weak shear seen towards the main contin- 
uum emission, a small scale shear is observed in the south- 
west edge of the field coinciding with weak 3mm continuum 
emission, that could there lead to the formation of a new 
protostar. 

We have therefore three MDCs (CygX-N12, N53, N48) 
which could correspond to three cores of similar nature 
seen at three different evolutionary stages: CygX-N48 cor- 
responding to a younger and CygX-N53 to a more evolved 
phase. CygX-N3 seems to be forming a cluster of at least 6 
protostars with the total mass shared between them with- 
out a dominating massive one. Finally CygX-N63 can be 
seen as a special case: a single collapsing object not dis- 
playing strong dynamics in a collapsing, large scale en- 
velope. Alternatively, it can be seen as the extreme case 
of the massive star forming MDCs discussed here (CygX- 
N12, N53, and N48) with only one protostar dominating 
the collapse. We therefore propose that the initial condi- 
tions may explain the differences between CygX-N3, CygX- 
N12/N48/N53, and CygX-N63 while the evolutionary state 
would be the explanation for the observed differences be- 
tween CygX-N12, N48, and N53. 

6. Summary & Conclusions 

We presented a detailed study of high-density gas tracers 
towards 5 MDCs in Cygnus-X focusing on their kinematic 
properties on global (~0.1 pc) and small scales (~0.03 pc). 
Our main findings are: 



1. On ~0.1 pc-scales dense gas traced by H 13 CO + and 
HCO+ shows infall signature in CygX-N12, -N48, -N53 
and -N63. The velocity fields indicate coherent, orga- 
nized bulk motions towards CygX-N3, -N53, -N63 with 
a rotation-like pattern, while towards CygX-N12 and 
-N48 no global velocity gradients are dominant. 

2. Simultaneous radiative transfer modeling of H 13 CO + 
and HCO + shows that due to depletion effects the emis- 
sion is tracing more the bulk of material around the 
protostars (~ 10 5 -10 6 cm -3 ). Depletion is found to be 
most significant in CygX-N63. 

3. At high angular-resolution we find a better correspon- 
dence of H 13 CN with continuum peaks than H 13 CO + . 

4. High angular-resolution maps with the PdBI in H 13 CO + 
reveal significant substructure at small-scales and in all 
MDCs several line components are found, which are dis- 
entangled into small-scale coherent flows with intrinsic 
velocity gradients and shears. 

5. Analysis of line-dispersion of the high- angular resolu- 
tion H 13 CO + spectra question if the level of micro- 
turbulence at small-scales is enough to provide sufficient 
support against gravitational collapse. 

6. At high angular-resolution the larger scale global coher- 
ent motions are resolved into individual coherent flows 
with filamentary structures, showing velocity shifts and 
more importantly shears/discontinuity in velocity to- 
wards the continuum peaks. The relative difference in 
velocity position of the flows give dynamical time-scales 
of the order of the free-fall time-scale for CygX-N3, - 
N12, -N48, -N53. Only for CygX-N63 we find an order 
of magnitude larger dynamical time-scales. 

7. Small-scale velocity shears may indicate ~10 4 yr dy- 
namical time-scales for new protostellar seeds built up 
by small-scale converging flows. 

8. Comparing the fragmentation and the kinematic prop- 
erties of our sample, they point to different evolutionary 
stages for CygX-N48, -N12 and -N53, while different 
initial conditions may explain the differences towards 
CygX-N3 and -N63. 
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Appendix A: Details on radiative transfer modeling 

To model the molecular line emission of MDCs, we consider 
a core embedded in its parent molecular clump with smooth 
transition of the physical parameters. For each core we used 
a simple model with the same assumptions: t he mass and 
the radius is taken from continuum data of Mot te et al.l 
( |2QQ7| ) and the corresponding 90% mass and radius were 
used to constrain the density profile with a power-law expo - 
nent of ai~-2 (similarl y as done in|Schneider et al.||2010 ). 
NH3 measurements by 



Wienen 



ml p008) give constrains on 
the temperature and the size of the NH 3 cores which cor- 
respond to the measured temperature. Using these param- 
eters as mass-averaged temperature within a size slightly 
larger than the size of the MDCs (~ 0.2 pc) the tempera- 
ture profile is then constrained. The cores are embedded in 
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a larger structure, a molecula r clump (^0.3-0.6 p c), whose 
parameters were taken from |Motte et aL ( 2QQ7[ ). We use 
the size of the clump for the outer radius and constrain the 
density profile with a shallower exponent (0.2 ~ -1.0 to - 
1.6). First estimates for turbulent widths are tak en fro m the 
gaussian line- fitting of the H 13 CO + lines. Figure A.l shows 
the density, temperature and abundance profiles used for 
the models. 

We use an abundance of X(HCO + )~10 -9 and intro- 
duce depletion at a given density. We both vary the den- 
sity at which depletion occurs (10 5 -10 6 cm" 3 ) and the level 
of depletion (0-1000). Summarizing, we keep the depletion 
density and the depletion factor as a free parameter and 
constrain the best fit using all the sp ectra e xtra cted along 



the radius (see the results in Figures A. 3 - A.7). Thus by 



fitting a radial profile, we restrict the degeneracy of the 
model. The explored parameter-space is the same for each 
MDC and the best fit is determined by comparing the re- 
duced x 2 values of the fit. From this modeling we obtain 
the opacity (r) and the density, where depletion occurs. 

Figure |A.2| shows a sketch of the model, the core em- 
bedded in a clump and the models of the spectra we get at 
certain radius. This shows the emission coming only from 
the core (a), then within the radius, where NH 3 observa- 
tions were taken (b) and finally the spectra taken through 
the whole clump (c) illustrating that a significant part of 
the emission is originating from 10 4 -10 6 cm -3 densities. 
Figure |A.1| shows the obtained density, temperature and 
abundance profiles, respectively. 





Log R [pc] 



Log R [pc] 



CygX-N1 2 
CygX-N48 
CygX-N53 
CygX-N63 



Log R [pc] 



Fig. A.l. Profiles of density, temperature and HCO + abun- 
dance. Colors indicate the different MDCs. 



CygX-N3 The HCO + spectra show red-shifted asymmetry 
in the line profiles, which may be an indication of expand- 
ing gas. Broad line- wings are present due to outflowing gas, 
which we do not aim to reproduce in our model. We con- 
firm with Simline that HCO + is optically thick and we 
obtain a reasonably good fit with an expansion velocity of 
~0.3 km s _1 to the central spectra. The radial fitting of the 
profile is less good indicating that we either overestimate 
the optical thickness of the line, or the expansion velocity 
is not constant. The spectral decomposition shows a sec- 
ond line component at the same velocity of the red-shifted 
peak, which may indicate, that expansion is not necessarily 
the reason for the line profile, but an optically thin second 




Fig. A. 2. The scheme of the physical model. Model spec- 
tra are produced for CygX-N3 in 3 variations, a) modeling 
only the core, b) the core corresponding to the NH 3 mea- 
surements, c) modeling the core and its embedding clump. 



component at the same velocity may produce similar line 
profile as well. Nevertheless it is beyond the scope of this 
work to perform such a detailed modeling. 

CygX-N12 Indications for contribution of more than one 
line-components can be seen in the spectra of CygX-N12, in 
addition to broad line- wings, which limits our approach to 
perform a good modeling for this core. Nevertheless our re- 
sults reproduce the line ratio of HCO + and H 13 CO + . Note, 
that the spectra show a shallower shoulder on the red part, 
which can be indication of infalling material. Therefore we 
include 0.2 km s _1 as infall velocity in the model. 



CygX-N48 and -N53 The clump associated with CygX-N48 
and C ygX-N53 was modeled in detail in Schnei der et al.| 
d20T0b. Here we use the same physical model for all the 
MDCs, which is adapted to model more precisely the cores 
e mbedded within the clu mp, that has not been considered 

polob. 



Schneider et al. 



CygX-N63 CygX-N63 has the weakest HCO+ and 
H 13 CO + emission compared to the other MDCs and the 
morphology of the high-resolution emission confirms a hole 
centered on the continuum peak. Our best model requires a 
high depletion factor and a lower depletion density (~2 10 5 
cm -3 ). Note that the line-profile is similar to CygX-N12, 
a shallow red shoulder of the spectra may indicate slight 
infall-motions that we model with 0.1 km s -1 infall 
velocity. 



Appendix B: Decomposition of spectral line 
components in H 13 CO + 

To disentangle the individual components seen in the spec- 
tral line maps (Figure [2] b) a double gaussian fit was per- 
formed to each spectra with an IDL routine using MPFIT. 
Figure [RT] shows the resulting two component fits in green 
and red. We recall that the noise towards the edge of the 
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position = 0" (0.000 pc) 



position = 14"(0.1 19 pc) 



position = 28"(0.238 pc) 



position = 42"(0.357 pc) 




Fig. A. 3. HCO + and H 13 CO + mean spectra computed over concentric annuli in steps of 14" with respect to the map 
center towards CygX-N3. The black line shows the observations, the red line corresponds to the model. 



position = 0" (0.000 pc; 



position = 14"(0.1 19 pc) 



position = 28"(0.238 pc) 




position = 42"(0.357 pc) 



Observation 
Model 




Fig. A.4. Same as Figure A.3 for CygX-N12. 



position = 0" (0.000 pc) 



position = 14"(0.119 pc) 



position = 28"(0.238 pc) 



position = 42"(0.357 pc) 




Observation 
Model 




Fig. A.5. Same as Figure A.3 for CygX-N48 



map is larger. Individual components are fitted in a fixed 
velocity range with respect to the vi sr . A line is fitted when 
the intensity is above S/N ratio of >5, otherwise no second 
component is fitted. It is a tendency that for CygX-N3, - 
N53 and -N63 the components are best separated in the 
outer parts of the core, while in the central part only a sin- 
gle gaussian is fitted. Especially, in core CygX-N63, most 
of the spectra are dominated by one strong component. 
This tendency is also consistent with the assumption that 
the double-peaked line-profiles are not due to an opacity 
effect (thus to self-absorption), because one would expect 
this to be significant in the positions of the strongest emis- 
sion associated with the central parts, which is clearly not 
the case. For CygX-N12 and -N48 the routine picks up well 
the two components across the whole field. Note that the 
line-positions of CygX-N48 are well corresponding to the 
position of the two components seen in H 13 CN (within ± 
0.3 km s -1 ). 



Table B.l. Result of spectral decomposition of IRAM 
30m+PdBI H 13 CO+ lines. 



Source 


peak velocity (vj sr ) 


line-dispersion (a) 




1st comp. 


2nd comp. 


1st comp. 


2nd comp. 




[km s -1 ] 


[km s -1 ] 


[km s -1 ] 


[km s -1 ] 


CygX-N3 


15.33 


16.75 


0.56 


0.58 


CygX-N12 


15.36 


16.87 


0.6 


0.68 


CygX-N48 


-4.46 


-2.75 


0.65 


0.65 


CygX-N53 


-4.43 


-3.17 


0.55 


0.54 


CygX-N63 


-4.39 


-2.89 


0.49 


0.36 



We summarized the fit parameters in Table |B.1| The 
identified components in the spectra correspond to an av- 
erage (where the line intensity is larger, than 5 S/N), since 
they are not fixed in position but show systematic shifts in 
velocity. 
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position = 0" (0.000 pc) position = 14"(0.119 pc) position = 28"(0.238 pc) position = 42"(0.357 pc) 



[ CygX-N53 j 














-10 -5 

km / s 



Fig. A. 6. Same as Figure A.3 for CygX-N53. 



position = 0" (0.000 pc) 



position = 14"(0.1 19 pc) 



position = 28"(0.238 pc) 



position = 42"(0.357 pc) 




Observation 
Model 



Fig. A.7. Same as Figure A.3 for CygX-N63 



Appendix C: Coherent flows seen at high spatial 
resolution - continued 

CygX-N12 shows several velocity components in the spectra 
and as mentioned in Section 5.2.4 due to their complexity 
it is not possible to determine their spatial structure in our 
systematic approach. Position- velocity cuts indicate a com- 
plex mixture of several line components and the extracted 
spectral profile confirm this complexity. Individual velocity 
components appear in a range of ~2.5-3.5 km s _1 . 

CygX-N48 shows a similarly complex mixture of indi- 
vidual velocity components with the individual components 
separated in a range of ~3.5 km s _1 . 

CygX-N53 may have a similar configuration as CygX- 
N3, but the separation of the individual components and 
thus the location of velocity shears is less clear. 
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Fig. B.l. Spectral decomposition of the line-components seen in 30m+PdBI H 13 CO + spectra maps. Black line shows tl§ 
observed spectra, while red and green lines show the resulting fits of two components. We only show the components 
with an S/N >5. 
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Fig. C.l. CygX-N12 presented similarly as Figure^ Integration range for the grey scale is 15.1 to 16.1 km s _1 , red 
contours show an integration range between 16.1 to 17.1 km s _1 , blue contours are integrated between 14.1 to 15.1 km s _1 . 
Contour levels go from 5xrms noise in steps of 3xrms noise. 
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Fig. C.2. CygX-N48 presented similarly as Figure [s] Integration range for the grey scale is —4 to —3 km s _1 , red contours 
show an integration range between —3 to —1 km s _1 , blue contours are integrated between —6 to —4 km s _1 . Contour 
levels go from lOxrms noise in steps of 3xrms noise. 
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Fig. C.3. CygX-N53 presented similarly as Figure [sj Integration range for the grey scale is —4.7 to —3.7 km s _1 , 
red contours show an integration range between —3.7 to —2.7 km s _1 , blue contours are integrated between —5.7 to 
—4.7 km s _1 . Contour levels go from lOxrms noise in steps of 3xrms noise. 
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